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Abstract 

We present multiply subtractive Kramers-Kronig (MSKK) rela- 
tions for the moments of arbitrary order harmonic generation sus- 
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ceptibility. Using experimental data on third-harmonic wave from 
polysilane, we show that singly subtractive Kramers-Kronig (SSKK) 
relations provide better accuracy of data inversion than the conven- 
tional Kramers-Kronig (K-K) relations. The fundamental reason is 
that SSKK and MSKK relations have strictly faster asymptotic de- 
creasing integrands than the conventional K-K relations. Therefore 
SSKK and MSKK relations can provide a reliable optical data inver- 
sion procedure based on the use of measured data only. 

1 Introduction 

In linear optical spectroscopy Kramers-Kronig (K-K) dispersion relations [?] 
have widely been exploited in data inversion both in absorption and reflec- 
tion spectroscopy [?] . The characteristic integral structure of K-K relations 
require the knowledge of the spectrum at a semi-infinite angular frequency 
range. Unfortunately, in practical spectroscopy only finite spectral range can 
be measured. Moreover, technical difficulties in gathering information espe- 
cially in nonlinear optical spectroscopy on a sufficiently wide spectral range 
make hard even the application of the K-K relation. The extrapolations in 
K-K analysis, such as the estimation of the data beyond the measured spec- 
tral range can be a serious source of errors [?, ?]. Recently, King [?] presented 
an efficient numerical approach to the evaluation of the K-K relations. Nev- 
ertheless, the problem of data fitting is always present in regions outside the 
measured range. 
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In the context of linear optical spectroscopy singly [?] (SSKK) and mul- 
tiply [?] (MSKK) subtractive K-K relations have been utilized in order to 
avoid the limitations due to finite-range data. Only recently they have been 
proposed in the context of nonlinear optics and especially for harmonic gen- 
eration susceptibility [?]. The idea behind the subtractive K-K is that the 
inversion of the real (imaginary) part of nth order harmonic generation non- 
linear susceptibility can be greatly improved if we have one or more anchor 
points, i.e. a single or multiple measurement of the imaginary (real) part for 
a set of wavelengths. Recently, we presented generalized K-K relations and 
showed their validity on polymer [?]. In this paper we present for the first 
time the theory for multiply subtractive generalized K-K for the moments 
of nth order harmonic generation susceptibility. In addition, we show using 
data of polysilane that the success of data inversion of the moments relations 
of third-order harmonic generation susceptibility is better already in the case 
of SSKK than for conventional K-K. The results of the latter were presented 
in Ref. [?]. This SSKK is important in the sense that usually it is rather 
difficult to obtain reliable information on nonlinear susceptibility at multiple 
anchor points. 
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2 Multiply subtractive Kramers-Kronig rela- 
tions for the moments of the harmonic gen- 
eration susceptibility 

We first present the theory related to moments of nth order harmonic gen- 
eration susceptibility. Then we deal with a function u> 2a x^{nuj; • • • s k>), 
where u> is the angular frequency, < a < n, and x^ n \ nuJ ) is t ne tensor 
element of the nonlinear susceptibility. From now on we make use of the 
abbreviation x^ n K nuJ 'i oj, - ■ ■ ,u>) = x { nu} )- Next we generalize the result of 
Palmer et al. [?] and Lucarini et al. [?] to hold also for the moments of the 
harmonic generation susceptibility. With the aid of mathematical induction 
(for a rigorous proof see Appendix A in Ref. [?]) we can derive multiply 
subtractive K-K relations for the real and imaginary parts of as follows: 




Re{x {n) (^i)} + "-- 



Re{x (n) (^)} + --" 




+ 2 [(^_^)(^_^)...(^_^)] P / 



oo 



,/2a+l 



Im{x {n) (nu')} 



Jo 



— u 



2 ) • • • K - <4) 



(1) 



4 




lm{ X {n) (nuj Q )} 



2a-l 



Im{x (n) (^i)} + --- 



lm{ X {n) (ncuj)} + ■ ■ ■ 



oo 



u' 2a Re{x (n) (nu')} 




with < a < n and P denoting the principal value integration. Here, ujj 
with j — 1, ■ ■ • , Q denote the anchor points. Note that the anchor points in 
Eqs. (UJ) and (0) need not to be the same. We observe that the integrands 
in MSKK relations presented in Eqs. (QJ) & (0) have remarkably stronger 
asymptotic decrease, as a function of angular frequency, than the conven- 
tional K-K relations. The integrands of the conventional K-K relations given 
for the Tith order harmonic generation susceptibility converge proportional 
to u~ ( - 2n+2 ^ at high angular frequencies, whereas the integrands in MSKK 
relations converge in relation to uj~ ( - 2n+2+2( ^\ Therefore, one can expect that 
the limitations related to the presence of an experimentally unavoidable fi- 
nite angular frequency range are thus relaxed, and the precision of the inte- 
gral inversion is then enhanced. According to Palmer et al. [?] the anchor 
points should be chosen inside the measured spectral range. Furthermore, 
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their study showed that accurate data inversion is possible when the anchor 
points are chosen near to the zeros of the Qth order Chebyshev polynomial 
of the first kind. In linear optical spectroscopy it is usually rather easy to get 
information of the optical constants such as real refractive index and absorp- 
tion coefficient of the medium at various anchor points. However, in the field 
of nonlinear optics it may be problematic to obtain the real and imaginary 
parts of the nonlinear susceptibility at various anchor points. Therefore, in 
the present study we wish to emphasize that the information even at a single 
anchor point reduces the errors caused by finite data range by utilization of 
SSKK. Now the choice of the location of the anchor point is not so critical. 
Nevertheless, the Chebyshev zeros accumulate at the ends of the spectral 
range. Therefore, we have chosen the anchor point near one end of the data 
range. From Eqs. (0) & (J2J) we obtain for the moments of the third-order 
harmonic generation susceptibility the following SSKK relations 

(3) 



7T J (u' 2 - U 2 )(u' 2 - uf) 

(4) 

where oj\ is the anchor point. 
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3 Application of SSKK relations to experi- 
mental data of third harmonic generation 
on polysilane 

We consider here the experimental data of Kishida et al. [?] related to third 
harmonic susceptibility of polysilane at energy range 0.4-2.5 eV. In Figs, la 
and 2a are shown the inverted data obtained using Eqs. Q & ©. In the 
case of Fig. la the anchor point was selected at the lower boundary of the 
data range. The data shown in Fig. la for the real part of the nonlinear 
susceptibility has an almost perfect match with the measured spectrum for 
< a < 2, while for a = 3 the finite range causes disagreement between the 
measured and retrieved data at low energy range of the spectrum. Neverthe- 
less, the spectrum is well retrieved at high energy range. In the case a = 4 
the theory does not predict convergence, but we have still reasonable agree- 
ment between experimental and inverted curve at the high energy range. If 
we compare the curves of Fig. la with those of Fig. lb obtained using the 
conventional K-K relations we observe that in Fig. lb the curve already for 
a = 2 is dramatically departing from the true one at the low energy range. 
Thus, the SSKK relations provide a better means than the conventional K-K 
relations to invert the moments data. 

In the case of the imaginary part of the nonlinear susceptibility the anchor 
point was selected at a location of the low energy side of the main spectral 
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feature, just out of the resonance. Also in the case of the imaginary part 
we observe that the SSKK procedure provides more precise data inversion 
than the conventional K-K relation (for comparison see Fig. 2b). In Fig. 
2a the agreement between measured and retrieved data is excellent for < 
a < 2 while for a = 3 we have a good agreement except for photon energy 
< 0.7 eV. In the theoretically non- converging a = 4 case, we somewhat 
surprisingly still have a good performance of the SSKK procedure at high 
energy range. Heuristically we can interpret this result as a manifest of the 
faster asymptotic decrease of the integrands realized when SSKK relations 
are utilized instead of K-K relations. 

4 Conclusions 

We have given, as far as we know, for the first time multiply subtractive 
Kramers-Kronig relations for the moments of the arbitrary order harmonic 
generation susceptibility. According to the calculations, which are based 
on experimental data of polysilane, we have observed how an independent 
measurement of the unknown part of the complex third-order (arbitrary- 
order) harmonic generation susceptibility for a given frequency relaxes the 
limitations imposed by the finiteness of the measured spectral data. The 
fundamental reason is that SSKK and MSKK relations have strictly faster 
asymptotic decreasing integrands than the conventional K-K relations. Thus 
SSKK and MSKK relations can provide a reliable data inversion procedure 
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based on the use of measured data only. That is to say no extrapolations of 
the data are needed beyond the measured range. We have demonstrated that 
even SSKK relations yield more precise data inversion, using a single anchor 
point, than conventional K-K relations. Naturally it is possible to exploit 
MSKK procedure if higher precision is required. However, the measurement 
of data at multiple anchor points may be experimentally tedious. 
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Figure captions 



Figure 1: Efficacy of SSKK relation (a) versus the conventional K-K relation 
(b) in retrieving Re{x^(3u;; uj, oj)} on polysilane. 

Figure 2: Efficacy of SSKK relation (a) versus the conventional K-K relation 
(b) in retrieving Im{%^(3(j;a;,a;,a;)} on polysilane. 
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Figure 1: Lucarini, Saarinen, and Peiponen 
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Figure 2: Lucarini, Saarinen, and Peiponen 
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